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The quasiparticle dispersion in the one-hole t — t' — t" — J model is studied. Both the finite-size 
diagonalization and the self-consistent Born approximation calculations have been performed and 
compared. The quasiparticle band structures in the hole and electron doped high-T c cuprates are 
qualitatively different. In the hole doped compounds, the band maxima locate at (±tt/2, ±7r/2), 
while in the electron doped compounds the band maxima locate at (n, 0) and its equivalent points. 
The angle-resolved photoemission data for the quasiparticle dispersion of Sr2Cu02Cl2 can be quan- 
titatively reproduced using the one band t — t' — t" — J model with the three-site hopping term. 



The energy dispersion of hole quasiparticles in the nor- 
mal states of two-dimensional CuO planes is of funda- 
mental interest for understanding the microscropic mech- 
anism of high temperature superconductivity. Due to the 
strong Coulomb repulsion between electrons, quasipar- 
ticles in high-T c cuprates behave very differently from 
what predicted by one-electron band calculations. Re- 
cently Wells et al. Q has reported an ARPES measure- 
ment on an insulating layered copper oxide S^CuO^C^. 
This is by far the most direct measurement of the disper- 
sion of a single hole in an antiferromagnetic background 
since Sr 2 Cu02Cl2 is difficult to dope. It provides a di- 
rect test for the model Hamiltonians which are proposed 
for high temperature superconductors. Wells et al. has 
found that the the hole bandwidth and the dispersion 
along the diagonal direction from (0, 0) to (tt, tt) in Bril- 
lion zone agree well with the calculations based on the t-J 
model. However, near (tt, 0) the data differ significantly 
from the prediction of the t-J model and the dispersion 
from (tt, 0) to (0, tt) is much greater than that in the t-J 
model. 

Recently Nazarenko et al. Q has calculated the hole 
dispersion in the one-band t — t' — J model under the 
self-consistent Born approximation (SCBA). They found 
that the presence of the i'-term can enlarge the energy 
dispersion around (0, tt), but the overall bandstructure is 
inconsistent with the experimental one. Several groups 
have also calculated the energy dispersions of hole in the 
three bands model based on either variational wavefunc- 
tions |U or SCBA @ and found that the results re- 
produce very well the experimental data. These studies 
seem to imply that the one band model is inferior to the 
multiband model even for studying low energy proper- 
ties of high-r c cuprates ||. This is actually not true. 
To understand qualitatively the physics of a hole prop- 
agating in an antiferromagnetic backgroun, the t — J or 
t — t' — J model might be sufficient. However, to fit 
quantitatively the result of a one band model with the 
experimental data, the one band model used should in- 
clude all the terms which are deduced from the multi- 
band Hubbard model by eliminating highlying orbitals. 
At least two of those terms which are ignored in the pre- 
vious study should be considered, one is the i"-term (i.e. 



the next-next neighbor hopping term) and the other is 
the three-site hopping term. Both the i'-term and t"- 
term are originated from the effective hopping (or wave- 
function overlap) between two O 2p orbitals besides a 
Cu via the Cu 3d and 4s orbitals ||; the i'-term is the 
hopping between two nearest neighboring Oxygens, and 
the £"-term is the hopping between two Oxygens on the 
two sides of Cu. \t"\ is generally smaller than but cer- 
tainly of the same order of magnitude as \t'\ |J. The 
three site hopping term is always present if the effective 
one band model is derived from the one band Hubbard 
model |?J or the multiband Hubbard model |8| . Both the 
t" and three-site hopping terms involve hoppings on same 
sublattices and are affected very weakly by antiferromag- 
netic correlations. They may therefore have substantial 
contributions to the quasiparticle dispersion. Other long- 
range hopping terms involve the wavefunction overlap of 
O orbitals not within the same unit cell and are generally 
small and negligible. 

In this paper we report our theoretical results for the 
energy dispersion in the one-band model with one hole in 
two dimensions. Both the Lanczos diagonalization and 
the SCBA calculations are performed and compared. We 
find that the experimental data for Sr2Cu02Cl2 can be 
quantitatively fitted by the t — t' — t"— J model with the 
three-site hopping term within the experimental error. 

Let us first consider the t — t' — t" — J model without 
the three-site hopping term. The model Hamiltonian is 
defined in the Hilbert subspace without double occupied 
sites by 
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where ( ) refers nearest neighbors, and tg = t for the 
nearest neighbor hopping (S — x, y), t' for the next near- 
est neighbor hopping (8 — x±y), and t" for the next-next 
nearest neighbor hopping (S = 2£, 2y). The rest of no- 
tation is standard. For high-T c materials, both t and 
J vary in a certain energy scale, t « 0.3 ~ 0.4ev and 
J sa 0.1 ~ 0.2ev, depending on compounds. So far ac- 
curate determinations for t and J from experiment are 
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still not available. For concreteness in discussion, we take 
t = 0.35ev and J = 0.15ev. The quasiparticle bandwidth 
evaluated from this set of values of t and J is roughly 
the same as that for Sr2Cu02Cl2 when t' = t" = 0. In 
our calculation, we take if and t" as free parameters, but 
restrict \t"\ < \t'\ < J. If the Hamiltonian (0) is de- 
rived from the multiband Hubbard model including the 
Cu 3d x 2_ y 2 and 4s and O p x and p y orbitals, then t' 
is negative (positive) for hole (electron) doped materials 
and t" has an opposite sign to if ||. 

For a one-hole system, we can do a Galilean transfor- 
mation to shift the origin of the frame of coordinates to 
the position of hole If the hole momentum is k, fol- 
lowing the derivation of Ref. || it can then be shown that 
the Hamiltonian (|l|) is equivalent to an effective Hamil- 
tonian 
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where = (0, 0) and P — J2ka ^ c \a Cka IS ^ ne total mo- 
mentum operator. In our exact diagonalization study, 
we solve this effective Hamiltonian on finite size lattices 
with periodic boundary condition (PBC), but allow k to 
change continuously (i.e. k can take all values allowed in 
an infinite square lattice). 

For each given fc, wc define the difference between the 
lowest energy of the one-hole Hamiltonian (Q) with total 
spin 1/2, EN-i(k), and the ground state energy of the 
system without holes, En, as the coherent hole quasi- 
particle energy E(k) = En — EN-i(k). The quasipar- 
ticle energy E(k) such defined corresponds directly to 
what measured in experiments, since the ground state 
of Sr2Cu02Cl2 is a spin singlet and the state with one 
electron removed from the ground state by high energy 
photons has spin 1/2. On finite size lattices, some higher 
spin states may have lower energies than EN-i(k) for 
certain values of k, so -E/v_i(fc) may not always be the 
minimum energy of Hh(k). However, in thermodynamic 
limit we believe that EN—\{k) will either be the min- 
imum of Hh{k) or differ infinitesmally small from the 
minimum energy of Hh(k) for all k. Thus alternatively 
E{k) can also be defined as the difference between the 
lowest eigeneigen of Hh(k) and the ground state energy 
without holes. In thermodynamic limit the above two 
definitions should give the same result for E(k). 

Under the SCBA |l(],[Il|], the quasiparticle energy is 
given by the position of coherence peak at the bottom of 
the spectral function A(k,uj) = —ImG(k,ui)/Tr, and the 
single particle Green's function G(k,u>) is determined by 
the self-consistent equation: 



G(k,u) = 
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where = At' cos k x cos k y + 2£"(cos 2k x + cos 2k v ) and 
T(k, q) = 4t(7 fe _ g u g + 7 fe w 9 ) with tt q = 4J y/1 - j q , u k = 



cosk y )/2. From Eq. (g), 



it is straightforward to show that G(k,u>) (and therefore 
the quasiparticle dispersion) is symmetric under the re- 
flection about the magnetic Brilliou zone boundary (i.e. 
k x + k y = 7i") since e k is symmetric under this transforma- 
tion. This symmetry is purely due to the approximations 
made in this approach; the original Hamiltonian (^) does 
not have this symmetry 

There are several approximations involved in the 
SCBA calculation: the linear spin-wave expansion and 
the neglect of crossing diagrams. These approximations 
ignore vertex corrections and the hole distortion to the 
spin background and do not guarantee the hard-core na- 
ture of the slave fermions and the Schwinger bosons. In 
the absence of the t' and t" terms, the contribution from 
the two-loop crossing diagram to vertexs is exactly zero 
|l2f . Thus the the vertex correction in the SCBA calcula- 
tion for the t-J model is small. However, in the presence 
of these terms, the two-loop crossing diagrams are gen- 
erally nonzero. 

For the t-J model, the band structure of hole quasi- 
particles has been extensively studied by many groups 
. and the agreement between the finite size calcu- 
lation jl3 14 and the SCBA calculation p2[ is remark- 
ably good. The quasiparticle dispersion in this model 
shows many interesting features which are completely dif- 
ferent from that in the ordinary metal. Firstly, the band 
maximum locates at (±7r/2, ±7r/2); secondly, the effec- 
tive mass is very anisotropic, it has a larger value along 
the zone diagonal and a smaller value along the direction 
perpendicular; and thirdly, the band width scales with 
a certain power of the exchange energy J (~ J 2 / 3 when 
i<C J) rather than the hopping constant t. 

In the previous finite size calculations, the Hamilto- 
nian (Q) with PBC were generally used. As the number 
of k points allowed in a finite lattice with PBC is lim- 
ited, calculations on large lattices are generally required 
in order to obtain a complete picture of the quasipar- 
ticle band structure. However, if the Hamiltonian (j|) 
with PBC is diagonalized and the hole momentum k in 
(0) is taken as a free parameter (this is eqivalent to di- 
agonalizing the Hamiltonian (fil) with twisted boundary 
conditions as shown in Ref. M ), we find that even on 
a relatively small lattice, such as N = 20, a comprehen- 
sive picture of the quasiparticle band structure can be 
obtained. 

In Fig. la, our finite size diagonalization results for 
E(k) on N = 16 and 20 are shown and compared with the 
SCBA ones for the t — J model. The agreement between 
the finite lattice calculation and the SCBA calculation in 
the whole Brilliou zone is good. The finite size effect, as 
revealed by the difference between the N = 16 and N — 
20 results and the difference between the finite lattice 
result and that of SCBA is small. On finite lattices, the 
band maxima locate not exactly at (±7r/2, ±7r/2), but 
they tend to move towards these points as increases 
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FIG. 1. Comparison of the quasiparticle dispersion relation 
E(k) (unit: ev) obtained from the SCBA method on 24 x 24 
lattice (solid line) with that obtained from the finite size diag- 
onalization method on TV = 16 (triangle) and N — 20 (circle) 
lattices for the t — t — J model with t = 0.35ev, J = 0.15ev 
and (a) t' = 0, (b) t' = -0.08ev , and (c) t' = 0.08ev. 

The presence of the next nearest neighbor and the 
next-next neighbor hopping terms changes largely the en- 
ergy dispersion of quasiholes. Let us consider the contri- 
bution of the t' term first. When t' < 0, as shown in lb, 
the energy dispersion along the line (tt, 0) — (0, 7t) is en- 
larged. Around (tt,0), the dispersion is relatively small, 
but E(k) falls much below the band maximum. These 
features are qualitatively consistent with the experimen- 
tal results for Sr2Cu02Cl2 PJ. However, the overall dis- 
persion of E(k) (except at the vicinity of (tt, 7r)) is largely 
reduced by the t' term. If fitting the experimental data 
of Sr2CuC>2Cl2 using the t — t' — J model, we find that 
a rather large J value, about twice as large as what has 
been found in experiments, is needed. 

When t' > 0, as shown in Fig. lc, E(k) behaves quite 
differently from the case t' < 0: (1) In this case the band 
maxima locate at (tt, 0) and its equivalent points; (2) the 
bandwidth is enlarged compared with the t— J model ; (3) 
E(k) becomes more symmetric about the magnetic Bril- 
lion zone boundary and shows a large dispersion along 
the line (tt, tt) — (tt,0) — (0,0); (4) The agreement be- 



tween the SCBA and the finite size calculations is fairly 
good and the finite size effect is even smaller than that 
for the t— J model. However, as t' gets larger, there is no 
well defined coherent peak at the bottom of quasiparicle 
spectrum when k is in the vicinity of (0, 0) in the SCBA. 
So far no ARPES data are available for electron doped 
cuprates for comparison. If we believe that the asymme- 
try between the electron and hole doped high-T c mate- 
rials is mainly due to the next nearest neighbor hopping 
term, then the above features of E(k) should in principle 
be observable in the ARPES spectrum in electron doped 
materials. 

For the t— t' — J model with a negative t' term, the over- 
all agreement between the finite size calculation and the 
SCBA calculation is not as good as for the t' > cases. 
This disagreement is probably due to the linear spin- wave 
approximation in the SCBA approach which ignores the 
influence of holes to the quantum antiferromagnetic Neel 
state. A positive (negative) t' term can enhance (reduce) 
the antiferromagnetic correlations in the one band model 
ifrsfl . So in the presence of a negative t' term, the dis- 
tortion of a hole to the quantum antiferromagnetic Neel 
state is enlarged. In this case to obtain a more quanti- 
tative description for the quasiparticle dispersion, exact 
diagonalizations on larger lattices are needed. 
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FIG. 2. Comparison of the quasipar- 

ticle energy E(k) - E(tt/2,tt/2) (unit: ev) in Sr 2 Cu0 2 Cl2 
(circle) with the corresponding results for the t — t' — t" — J 
model obtained from the exact diagonalization with N = 20 
(solid line) and the SCBA approach (dashed line), t = 0.35ev, 
t' = -0.12ev, t" = 0.08ev, and J = 0.15ev. 

The above discussion indicates that the bandwidth of 
the t—t' — J model is too small compared with the experi- 
mental one for Sr2Cu02Cl2 within the physically reason- 
able region of parameters, consistent with the previous 
study ||. A finite t" term, however, can change this situ- 
ration significantly. The contribution of the t" term to 
E(k) is similar to a —t' term, for example for the t—t" — J 
model (t' = 0) the band maximum locates at (tt/2,tt/2) 
if t" > or (tt, 0) if t" < 0, but the overall energy disper- 
sion is always enlarged by the t" term. 

In Fig. (g) the quasiparticle dispersions for the t — t' — 
t" — J model are shown and compared with the experi- 
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mental data for Sr2Cu0 2 Cl 2 . Along two diagonal lines, 
(0,0) — (tt, 7r) and (tt, 0) — (0,7r), the finite size effect is 
small and the agreement between the experimental data 
and our calculation is very good. Compared with Fig. 
lb, we find that the bandwidth is largely enhanced by 
even a small i"-term. Along the line (tttt) — (tt, 0) — (0, 0), 
the energy dispersion obtained from both the exact di- 
agonalization and SCBA calculations is larger than the 
experimental one, and the agreement between the exact 
diagonalization and the SCBA results along this line is 
also not as good as along two diagonal lines. A sub-peak 
appears in the curve of the exact diagonalization result 
on the line (tt, 0) — (0, 0). This is purely due to the finite 
size effect. 




It implies that the the multiband Hubbard model 
is indeed equivalent to a one-band model in describing 
low energy excitations of high-T c cuprates. 

In conclusion, we have studied the energy dispersion 
of hole quasiparticles in an antiferromagnetically corre- 
lated background using the exact diagonalization and the 
SCBA methods. In the exact diagonalization study, we 
first derive an effective Hamiltonian for the one-hole t-J 
model using the Galilean transformation for each given 
hole momentum k (which can take all allowed values in 
an infinite lattice), and then diagonalize this effective 
Hamiltonian on finite lattices with PBC. As no limita- 
tion for k values, this allows us to have a comprehen- 
sive study for the band structure of quasiparticles with 
the finite size diagonalization even on a N — 20 lattice. 
We find that the finite size diagonalization results agree 
well with the SCBA ones, especially when t' > 0. For 
hole doped high-T c compounds, the band maxima locate 
at (±7r/2, ±7r/2) in thermodynamic limit. However, for 
electron doped compounds, the band maxima locate at 
(tt, 0) and its equivalent points. The experimental data 
for Sr 2 Cu02Cl2 can be quantitatively understood from 
the one band t — t' — t" — J model with the three-site 
hopping term. 
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FIG. 3. Comparison between the quasiparticle energy 
E(k)-E(w/2,w/2) (unit: ev) of the t — t' — t" — J model with 
the three-site hopping term on a N = 20 lattice (curve) and 
that of Sr 2 Cu0 2 Cl 2 . t = 0.35ev, t' = -0.12ev, t" = 0.08ev, 
and J = 0.15ev. 

Now let us consider the three-site hopping term. The 
Hamiltonian for the three-site hopping term is given by 
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This term describes an effective hopping of a hole to one 
of its next or next-next neighbor sites by exchanging spins 
with another hole on its nearest neighbor site. From the 
calculation we find that this term has very weak effect 
on the quasiparticle dispersion along two diagonal lines. 
However, it has a relatively larger effect on E(k) when 
k varies along (tt, tt) — (tt, 0) — (0,0). It suppresses the 
dispersion of E(k) on (tt, tt) — (tt, 0) — (0, 0) and therefore 
improves the agreement between the theoretical result 
and the experimental data. Fig. (|^) compares the energy 
dispersion of the t — t' — t" — J model with the three-site 
term on a N = 20 lattice with that of Sr2Cu02Cl2- It is 
easy to see that the one-band t — t' — t" — J model with 
the three-site hopping term gives a good account for the 
quasiparticle dispersion in Sr2Cu02Cl2. This result is 
consistent with the previous studies for the quasiparti- 
cle dispersion based upon the multiband Hubbard model 
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